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We review the analytic results for the phase shifts 6i(k) in non-relativistic scattering from a spherical well. The 
conditions for the existence of resonances are established in terms of time-delays. Resonances are shown to exist for 
p-waves (and higher angular momenta) but not for s-waves. These resonances occur when the potential is not quite 
f~^. ' strong enough to support a bound p-wave of zero energy. We then examine relativistic scattering by spherical wells and 
barriers in the Dirac equation. In contrast to the non-relativistic situation, s-waves are now seen to possess resonances 
in scattering from both wells and barriers. When s-wave resonances occur for scattering from a well, the potential is 
not quite strong enough to support a zero momentum s-wave solution at E — m. Resonances resulting from scattering 
from a barrier can be explained in terms of the 'crossing' theorem linking s-wave scattering from barriers to p-wave 
scattering from wells. A numerical procedure to extract phase shifts for general short range potentials is introduced 
. and illustrated by considering relativistic scattering from a Gaussian potential well and barrier. 
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Qh' The question of what does or does not constitute a resonance in potential scattering still seems to cause some 
confusion. Many authors who consider scattering in simple systems such as the one-dimensional square well in the 
Schrodinger equation talk of transmission 'resonances' when the transmission coefficient T = 1 with little or no 
justification as to why these might be resonances. It would seem on the surface that a resonance exists on account of 
a maximum or peak in a measure of the scattering (for example the transmission coefficient in one dimension or the 
cross-section for a particular partial wave in three dimensions) . But Wigner [1] pointed out nearly 50 years ago that in 
order to have a true physical resonance the scattered particle must be captured in some way by the scattering centre. 
So in addition to the peak in the scattering probability function, there must also be a delay in the transition time 
of the particle through the target: this requirement is fulfilled by demanding that the Wigner time delay is positive. 
This implies that the phase shift associated with the peak must increase through an odd multiple of ^ as the particle 
energy increases through the resonance energy. 

We will begin this paper by looking at the phase shifts for the simple case of scattering from a spherical well * in the 
Schrodinger equation in order to illustrate resonant behaviour and distinguish it from cases where the cross-section 
peaks but there is no resonance. The argument and results closely follow those given by Newton [2]. It will be seen 
that, using Wigner's criteria, the spherical well does not give rise to s-wave resonances. Resonances do exist, however, 
for scattered p-waves and waves of higher angular momenta. This comes about because of the centrifugal term in the 
Schrodinger equation for partial waves I > 0. 

The relationship between bound states and continuum states for non-relativistic systems is well documented. Schiff 
[3] noted that 'a potential well that has an energy level nearly at zero exhibits a resonance in the low-energy scattering 
of particles with the same I value as the energy level'. Newton [2] later studied this more carefully and showed that 
one can go smoothly from a zero energy solution to a real resonance by weakening the potential slightly. This idea 
is crucial to the the procedure we will adopt. First the condition on the potential to support a zero energy solution 



"this is often referred to as the three-dimensional square well 
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must be established. This condition can then be relaxed slightly to push the solution into a continuum state. Finally 
the phase shift corresponding to the resulting peak in cross-section is checked to see if it fulfils the criterion of positive 
time delay necessary for a real resonance. 

In this paper we generalise these results to the relativistic Dirac equation in three dimensions. We begin by 
reproducing the results for the phase shift in the spherical potential well and analyse whether resonant behaviour 
exists. In contrast to the non-relativistic results it will be shown that scattered s- waves do now give rise to resonances. 
It will be shown that these arise because of the presence of a centrifugal term in the original coupled equations which 
does not disappear even for s-waves. But we cannot stop there: in the Dirac equation we must consider the zero 
momentum solutions which exist at both E — m and E = —m. We refer to these as critical solutions and in particular 
it is conventional to refer to the solution at E = —m as supercritical. It will be shown that a 'crossing theorem' exists 
which connects scattering of particles by potential barriers to the supercritical state at E = —m. This phenomenon 
was first illustrated in a previous paper by two of us [4] in the context of positron-heavy nuclei scattering. We use this 
theorem to show that resonances exist in s-wave Dirac scattering from spherical barriers as well as spherical wells. 
Resonances also occur in higher partial waves as expected. 

We then investigate resonant scattering by general short range spherically-symmetric potential wells and barriers 
in the Dirac equation. A numerical procedure for the extraction of the phase shifts is demonstrated in the case of a 
Gaussian potential wells and barriers, and the results are analysed. We demonstrate that s-wave resonances can exist 
in scattering by both barriers and wells. 



THE SCHRODINGER EQUATION 

The phase shifts 8i(k) for the spherical well V(r) = —V r < a; V(r) = r > a are well known [3], [5]: 

tangi = kj[{ka)ii{pa) - p ji{ka)j[(pa) 
1 kn[{ka)ji(pa) - pni(ka)j[{pa) 

where ji and Hi are the regular and irregular spherical Bessel functions respectively, k 2 — 2mE and p 2 = 2m(E + V). 
So for s-wave scattering when I = this becomes: 

k tan pa — p tan ka 

tan d = ; ; (2) 

p + k tan ka tan pa 

while for p-wave scattering when I = 1 this gives: 

akp 2 tan pa — ak 2 p tan ka + (k 2 — p 2 ) tan ka tan pa 
tanoi = — — — ; — - : (3) 

ak z p + (p z — k z ) tan pa + akp z tan ka tan pa 

The potential wells of the greatest interest are those which are not quite deep enough to support a bound state. The 
critical values of the potential well which support states of zero energy (these may or may not be bound - this will be 
discussed later) are well established [3] . For ( = 0we have 

, , , , (2n - 1) 2 tt 2 
jo (pa) =paji{pa) => cos pa = =>V C = (4) 

and for I > we have 



For I = 1 this leads to 



ji-i(pa) = (5) 

2 2 

smpa = =► V c = — 2 (6) 

We are now in a position to examine the low energy phase shift behaviour in order to see if resonances occur. In the 
case of the s-wave, if the potential is not quite strong enough to bind the first state, the phase shift 5o(k) will be seen 
to rise from zero but never quite reach the value ^ before it decreases for larger energies. If we continue to increase 
the strength of the potential well there will be a value at which the phase shift at zero energy So(0) flips from zero 
to |. The strength at this point is exactly the critical value V — V c from Eq. (4), required for the potential to just 
support the first I = critical state. As soon as we just increase the strength of the potential from V c , <5o(0) flips to 
the value n and decreases for increasing energy. This behaviour is illustrated in the following figure: 
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FIG. 1. The s-wave phase shifts for the spherical well. A zero-energy state exists for a potential of strength V c = 0.257T 2 « 2.47 
(dotted line). The solid line represents a potential of strength V = 2.4 which almost binds whilst the dashed line is the potential 
of strength V = 2.6 which contains a loosely bound state, (a = 1 and m = 0.5). 



The process repeats with every additional bound state, so <$o(0) increases by ir each time and we have Levinson's 
theorem: 

<Jo(0) - n Q TT (7) 

where n tells us the number of I = bound states. For critical values of the potential when it is just strong enough 
to possess a zero-energy state, Levinson's theorem must be modified to 

S (0) = (n + ±)tt (8) 

Although it is not obvious, the critical state is not actually bound but instead forms part of the continuum. In order 
to see this consider the Schrodinger equation with for a potential well: 



tf"(r) 



fcil + 2mV(r) - k 2 



*(r) = (9) 



If the potential vanishes faster than 1/r 2 at large distances, the zero-energy solution for large r is found from 

*» « ^r^M ( 10 ) 

Consequently <J> behaves as 1/r' at infinity and the particle will only be normalisablc, hence bound, for I > 0. This 
means that for the first critical value of the potential, V = V c = 7r 2 /4 when I = 0, the state is not bound, no = and 
5o(0) = | from (8). These states are often referred to as half-bound states [2]. 

In order for a resonance to exist the Wigner time delay must be positive. This is defined as 

T ? dSl 2d6t mi 
Tl = 2 dE = vlk 

So the phase shift must increase through 7r/2. A decrease in the phase shift as the energy increases through ir/2 gives 
rise to a negative time delay - a time advancement [6]. This does not constitute a resonance. Hence there are no 
s-wave resonances for a spherical well. 

The results for p- waves (and other higher angular momentum waves) are distinctly different. If the potential is 
not quite strong enough to bind the first state, the phase shift 5i(k) rises from zero but it now increases through 
the value ^ at some positive value of the energy up to a maximum value (< n) before decreasing for larger energies. 
If we continue to increase the strength of the potential, the maximum of the phase shift approaches the value 7r for 
decreasing values of the energy until the phase shift at zero energy 5\(0) flips from zero to ir. The strength of the 
potential at this point is exactly the critical value V c = n 2 from (6), required for the potential to just support the first 
I = 1 bound state. It should be noted that, contrary to the s-wave case, the critical wave functions are bound states. 

To illustrate that there are p-wave resonances but no s-wave resonances, we can consider the effective potential 
V eff (r) = V(r) + 1(1 + l)/r 2 shown in figure 2. 
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FIG. 2. The effective potential for a spherical well. 



For I — there is no centrifugal barrier to trap the particle for small positive energies and consequently low 
momentum resonances do not occur. At zero energy there is nothing to stop the particle escaping the well - it is not 
bound and forms part of the continuum. For angular momenta I > 0, there is now a barrier which can momentarily 
trap particles giving rise to a resonance in the physical region E > 0. If we scatter a particle of energy E < 1(1+ l)/a 2 
(i.e. below the lip of the repulsive centrifugal barrier - see figure 2), it can tunnel into the target and once there 
eventually leak back out. If we deepen the potential the particle becomes increasingly trapped as the strength of 
the barrier through which it must tunnel increases and the resonance becomes more pronounced as it is retarded for 
greater lengths of time. We finally arrive at the critical value of the potential for which E = and the particle is 
completely trapped provided / > 0. 

As we continue to increase the strength of the potential beyond this critical value 5/(0) jumps by ir for each additional 
bound state as is the case for I = 0. For I > there are no half-bound states so Eq. (8) is no longer applicable and 
Eq. (7) can be modified to 

<5/(0)=n/7r (12) 

where n; now tells us the number of bound states with angular momentum I. From figure (3) we can see that the 
phase shift associated with potentials which almost contain a bound state will increase through an odd multiple of 
■| . As we get closer to criticality the phase shift increases more and more rapidly through this value and will peak at 
values that tend to (mod ir). This manifests itself as an increasingly more pronounced peak in the cross-section. 

This behaviour is illustrated in figure 3. The solid line represents a potential of strength V = 9 and the dashed 
line represents a potential of strength V = 9.7, where both potentials are just too weak to have a bound state. The 
dotted line is a potential of strength V = 10 which contains a loosely bound state. A zero-energy I = 1 bound state 
exists for a potential well of strength V c — ir 2 ~ 9.87. (a = 1 and m = 0.5). 
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FIG. 3. The spherical well for p-wave phase shifts. 



THE DIRAC EQUATION : SCATTERING SOLUTIONS 

We will now consider the problem of scattering from a spherical well in the relativistic Dirac equation. Following 
Grcincr, Miiller and Rafelski [7], which we refer to as GMR, the coupled radial equations for an electron of mass m 
and energy E in the presence the spherical well V(r) = —V r < a; V(r) = r > a can be written 
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f'(r) = -±f(r) + (E + m + V)g(r) 



g'(r) = -(E-m + V)f(r) + *g(r) 



(13) 



where the Dirac wave function *(r) for the radial equation is given by r*(r) = ^ ^ ^ The variable x = ±(j + 5) = 
±1, ±2, ... with the orbital angular momentum I = j + ={0, 1, 2, ..} corresponding to {s,p, d, ..} waves. When 



r < a, the solutions are 



/(r) = airji x (pr) 



where the well momentum p is 



\ X \E + V + m 
p 2 = {E + V) 2 - m 2 



31 -Aw) 



and 



X 

-x-i 



I x-i 



-x>o 
-x<o 



(14) 



(15) 



(16) 



When r > a, the free particle solutions are: 



/(r) = hrji x (kr) + b 2 rni x (kr) 
g(r) = — - ; [&iji_ x (fcr) + 6 2 n;_ x (fcr)] 



|X| S + m 

where k 2 = E 2 — m, 2 . When r — > 00 the following asymptotic forms for the spherical Bcsscl functions are required: 



(17) 



xj n {x) ~ sin (a; - ^f 1 ) 
The components can therefore be written as 

f{r) 

9(r) 



xn n (x) 



~ — cos x 



YJ 



61 sm kr ^— 

2 



1 

k 

X 1 

|X| S + m 
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6 2 cos [ fcr 



61 sin I kr — ^ 



— b 2 cos I kr — * ?r 



Therefore 



/(r) 
.9(0 



where 



So 



61 4 
]h\k 
b i X A 
~* |6i| \ X \E + m 

A cos 5, = bi 



tan <5; = - — 
61 



sin ( £;r - + Si 



sin I fcr — * n + Si. 



Asm. Si = -b 2 



A 2 = b\ + b\ 



(18) 

(19) 
(20) 

(21) 
(22) 

(23) 
(24) 



This relationship between the phase shifts and the coefficients 61 and b 2 was shown for s-waves by GMR. Depending 
on the sign of x we have 



_ / l x + 1 x < 
l x -\ x > 



X 



(25) 
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Using the above together with Eqs. (21 and 22), the two components are found to be 

bi A . f l l x ir 



\h\E + m V 2 
The asymptotic forms are clearly dependent on the sign of %: 



^Mr in r -2 + N (26) 

5(0^1^T7t4— cos (^-^ + (27) 



fM\ 6i / x cos (Ar + ^ + 
rtt(r) = ; r - * X < (28) 




f(r) \ h ( i**(kr-*£+Si x ) 



r*w = ( ;;;; )-^VH i ,. ^ . , I *>° (29) 



These results are consistent with Barthelemy [8], [9]. We can combine Eqs. (17) and (24) to incorporate the phase 
shift into the free-particle wave components of the wave function: 

f(r) = b 1 r[ji x (kr) - tan Si x ni x (kr)} 

9{r) = ^\ET^ bl[jl -* {kr) ~ t!m K n i-A kr )] ( 30 ) 

The solutions for r < a from Eqs. (13) and r > a from Eqs. (30) must be joined smoothly; by dividing the top 
component by the bottom component and equating the two solutions at r = a we arrive at 



aiji x (pa) bi \ji x (ka) - tan Si x m x (ha)] 



After re-arrangement this becomes 

k E + V + m ji x (pa) ji x (ka) - tan 5i x n lx (ka) 
p E + m ji_ x (pa) ji_ x (ka) - ta,nSi x ni_ x (ka) 



(31) 



(32) 



Taking p = ^(E + V) 2 — m 2 we find that the kinematic term 7 is given by 

P E + m 

1= k E + V + m (33) 
Solving for tan<5; x we find that the relativistic equivalent of the Schrodinger result (1) is 

tan5; _ 13i x (ka)ji_ x (pa) - ji x (pa)ji_ x (ka) 
lx 13l- x (P a ) n i x ( ka ) - 3l x (pa)ni_ x (ka) 



THE DIRAC EQUATION : ZERO MOMENTUM SOLUTIONS 

In this section we establish the critical solutions for the spherical well. Unlike the previous non-relativistic situation, 
the Dirac equation has two critical solutions corresponding to E = m and E = —to; the analysis which follows allows 
us to consider both of these critical solutions. For the spherical well potential the coupled radial Eqs. (13) have the 
solutions inside the well: 

f(r) = r[aiji x (pr) + a 2 m x (pr)} 

ffW = U\ E + V + to [aiJl - x {W) + ° 2n '-* (pr)] (35) 



G 



and the solutions outside the well arc 



/(r) = r 



2nr 



IhKi + i{nr) +b 2 Ii + i(Kr)} 



9{r) = eT^ \IW [ - blKl -^ M + {Kr)] 



(36) 



where k 2 



E 2 . K n+ i(nr) and I n+ i(nr) arc the modified spherical Bessel functions. In order to establish 



the bound states the internal wave function must be normalisable at the origin therefore a 2 = 0. The external wave 
function must also be normalisable at infinity which imposes the condition 62 = 0. By matching the top and bottom 
components at r — a the following conditions are obtained: 



2na, 



aiji x {pa) = \j—hK lx+ i (no) 



X 



P 



aiji_ x {pa) 



\ X \E + V + m 
Dividing the first equation by the second gives: 



Ixl P3i- X (pa) 



Upon re- arrangement this becomes 



E + m V 7T 



= -(E + m) 



2Ka h is 1 \ 



K, , 1 (na) 
nK t + i(na) 



3U (pa) 



X P E + m 



K, 



l +i 



(na) 



ji- x (P a ) \x\ k E + V + mK l x+ i(Ka) 
The critical limits E — ► ±m imply k — > 0, which together with the following identity from reference [10] 

~2 



(37) 



(38) 



(39) 



lim 



K n+h {z) 



{2n-l)\\z- n -\ 



(40) 



allow us to find the zero momentum limit of Eq. (39). Using Eq. (15) with E — ±m it is easy to see that 
p^p± = \/V(V ± 2m). Eq. (39) becomes 



3i x {P±a) 



X P± E + m (2Z X - 1)! 



ji_>±<*) Ixl « E + V + m(2l^ x -iy.\ 

We now consider x < and x > separately. 



(na) 



-{lx-l-x) 



(41) 



Case I : x < 



When x < this simplifies to 



.i-x-i(P± a ) P± a E + m 



j- x (p±a) 1 + 2 X E + V + m 

For critical states at E = m, the above can be written as 



E^±rn, k^0 



3-x-i(P+ a ) 2ma / V 



j- x {p+a) l + 2xVU + 2m 



(42) 



(43) 



The supercritical states at E = —to satisfy 



3-x-i(P- a ) = 



(44) 
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In particular for critical s-waves (x = — 1, l x = 0, l- x = 1) at E = m, Eq. (43) is 



ji(p+a) " V V + 2m 
Upon simplification this reduces to the transcendental equation 



tan a v /V r (^ + 2m) = -2m fl y - ^ 2 — (46) 

Supercritical s-waves satisfy jo(p_o) = sin(p_a) = from Eq. (44). This imposes the following condition on the 
well depth: 

a \/V{V - 2m) = nir (47) 



Case II : x > 

When x > 0, Eq. (40) simplifies to 

3 x (p±a) _ P±(2x-1) 1 



(48) 



j x -i(p±a) a (m-E)(E + V + m 

The critical solutions E = m give 

j x -i(p+a) = (49) 

while the supercritical solutions at E = —m are 



J X (P-°) 2 X -l,/y-2m (5Q) 



j x _i(p_a) 2ma V V 

For critical p-waves ( X = 1, Z x = 1, i_ x = 0) at E = m, Eq. (49) imposes the condition j (p + a) = => sin(p + a) = 0. 
The potential depth must therefore satisfy 

aV^(U + 2m) = wr (51) 

p-waves become supercritical when 



jo(p_a) 2ma V V 
Upon simplification this becomes 



jiip-a) _ 1 ,/U-2m ({ . 2) 



, V — 2m 

tanaVU(U - 2m) = 2tooW — (53) 

These conditions can also be derived from the bound state eigenvalue equations provided in GMR. 



PHASE SHIFT BEHAVIOUR : THE SPHERICAL WELL 



In this section we start by considering the conditions on the potential for a resonance to exist when low-momentum 
electrons are scattered on a spherical well. These conditions will then be shown to be related to the results for the 
critical solutions derived in the previous section. Finally the s-wave phase shift is explicitly derived and both s- and 
p-wave phase shifts are illustrated for spherical wells which nearly support an s or p zero momentum solution at 
E = m. 
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A necessary condition for the existence of a resonance is that the phase shift is Si x — ^. Further analysis is required 
to establish if this position corresponds to an increasing phase shift to give a resonance. In order to ensure this the 
denominator on the right hand side of Eq. (34) must be zero which gives the condition 



3l x (P a ) 



7: 



ni (ka) 



For low momentum scattering E 



jl_ x (pa) ' m_ x {ka) 
m on the spherical well p — > p+ and 
\im[x l+1 rn(x)} -> -(21- 1)!! 



Eq. (54) becomes 



Using Eq. (33) 



Eq. (56) becomes 



When x < 0, this is 



jijp+a) 
ji- x (p+a) 



— kaj 



~ 1)!! 
(2Z_ X - 1)!! 



(ka) 



kaj = ap + 



E- 



E + V + m 



2ma 



-Ox-'-x+i) 



V 

V + 2m 



2ma 



V 



(2^ 



1)!! 



K (p+ a ) 

jl_ x (p+a) ~'~~V V + 2m (2/_ x - i 



(ka) 



-Vx- 



2ma 



j-x(P+ a ) 



V 



1 + 2\ V V + 2m 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



which is identical to Eq. (43), the condition for x < critical states to be at E — m. When x > 0, the right hand 
side of Eq. (58) is infinite leading to the condition 



j x -i(p+a) = 



(60) 



which is identical to Eq. (49), the condition for \ > critical states to be at E — m. The implications of these results 
will now be highlighted for s- and p-wave scattering. 
For s-waves Eq. (34) simplifies to 



tani5 = 



ljo(ka)j 1 (pa) - jo(pa)j 1 (ka) 



jji(pa)n (ka) - j (pa)ni(ka) 
Using standard identities for ji(x) and ni(x) with I = 0, 1 we find: 

(jk — 1) tan /ca tan pa — jkpat&nka — fcatanpa 



tan <5 



akp tan ka tan pa + (1 — jk) tan pa + jkpa 



(61) 



(62) 



If we compare this to the non-relativistic case it can be seen that this more closely resembles the p-wave phase shift 
(Eq. 3) than the corresponding s-wave phase shift (Eq. 2). The difference in behaviour of this s-wave phase shift is 
clear in the figure below: 



So 
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FIG. 4. The s-wave phase shifts for a Dirac particle in a spherical well. The critical potential V c = 4.20m. The solid lines 
represent potentials of strength V = 4.195m (middle line) and V = Am (lowest line). The dashed line is a potential well of 
strength V = 4.3m which contains a bound state, (a = 1/m). 



It can be seen that the s-wave phase shift now just increases through | for values of the potential that are weaker 
than a critical value V c and thus we have a very weak resonance by our definition above. As will be shown below, 
the width of this resonance is large compared to that for the p-wave. This is a new result - we would not expect to 
obtain resonances for s-waves in a spherical potential well. 

The p-wave phase shifts are very similar to those obtained for p-waves in the Schrodinger case and are illustrated 
below. 




FIG. 5. The p-wave phase shifts for a Dirac particle in a spherical well. The critical potential V c = 2.30m. The solid lines 
represent potentials of strength V = 2m (lowest line) and V = 2.2m (middle line). The dashed line is a potential well of 
strength V = 2 Am which contains a bound state, (a — 1/m). 



For energies close to the resonant energy position, Er, the resonance has the Breit-Wigner form: 
where T; x is the width. This width is related to the time delay [11] by 



(63) 



r» . ~ 



(64) 



By taking the ratio of the widths of the s-wave resonance for V = 4.195m and the p-wave wave resonance for 
V = 2.25m, the s-wave resonance is seen to be approximately 230 times wider. 

For a very clear illustration of the existence of s-wave resonances for the relativistic spherical well, we use the values 
for the well width a and depth V of Pieper and Grcincr [12] in their discussion of superheavy nuclei. In that paper 
they establish the supercritical values for s states but do not consider scattering phase shifts. Here we choose a to 8 
fm and then the first three critical values at E = m and E = — m for s and p states are easily calculated using Eqs. 
(50, 51, 55 and 57) and are tabulated below: 



s (E=m) 


s (E=-m) 


p (E=m) 


p (E=-m) 


75.947 
153.434 
230.919 


77.997 
155.480 
232.964 


76.975 
154.458 
231.942 


79.012 
156.498 
233.983 



Table 1: The critical well depths in MeV for s- and p- states at E = ±m. 

The critical well for the Is i state is now weakened slightly by 1% to make the depth V = 75.187 MeV. The phase 
shift is seen to be: 
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E[keV] 



0.5 1 1.5 2 2.5 850 900 950 1000 

FIG. 6. The lsi phase shift and energy dependence of sin 2 So for the spherical well (a = 8 fm and V = 75.187 MeV). 



The phase shift is clearly seen to increase through -| corresponding to a resonance. The energy dependence of 
sin 2 So is plotted beside the phase shift and the width is found to be T = 41.83 keV. 

Similarly the critical well for the 2pi state is weakened slightly by 1% to make the depth V — 76.205 MeV. The 

phase shift and energy dependence of sin 2 Si is seen to be: 



Si 



k[m] 




E[keV] 



0.5 1 1.5 2 2.5 850 900 950 1000 

FIG. 7. The 2pi phase shift and energy dependence of sin 2 81 for the spherical well (a = 8 fm and V = 76.205 MeV). 



The width of the 2pi resonance is seen to be Ti = 11.5 keV. The 2pi phase shift is seen to increase more rapidly 
through t£ than the lsi phase shift and consequently the lsi resonance is wider than the 2pi resonance by a factor 
of approximately 3.6. 

In order to understand why s-wave resonances occur for relativistic particles we de-couple the large, /, and small, 
<?, components of Eqs. (13) to form two second order equations. In the case of a general short-range potential the 
equations for solutions outside the potential are: 



.9" - 



£2_ m 2_X(X + l) 



p2 2 X(X-l) 

E — m - 



g(r) = 
f(r) = 



(65) 



What happens is that the angular momentum barrier for the large component vanishes for si -waves when x = — 1 
but the small component still has an angular momentum term. In contrast for pi-waves (x = 1), the barrier remains 
for the large component but disappears for the small component so the centrifugal barrier never vanishes for both / 
and g. With this angular momentum term in one or more of the components wc arc considering cases analogous to 
the I > type scattering in the Schrodingcr equation and hence this gives rise to real resonances. 



PHASE SHIFT BEHAVIOUR : THE SPHERICAL BARRIER 



In the case of low momentum scattering E — > m on the spherical barrier, the analysis of the previous section can 



be repeated with V — > — V. The barrier momentum is now p —> q = y/{E — V) 2 - 
becomes q — > p_ . The corresponding form of Eq. (59) for scattering from barriers is 



m 2 . In the limit E — > m this 
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j-x-AP- a ) 



2ma 



V 



j_ x (p_a) l + 2x\V-2m 
which is the supercritical result for x > given by Eq. (50). Similarly Eq. (60) becomes 

j x -i(p-a) = 



(66) 



(67) 



which is the supercritical result for particles with x < given by Eq. (44). 

In the previous section we showed that the conditions for s-wave resonance scattering at low momenta off a spherical 
well agreed with the critical conditions at E = m calculated from the bound s state spectrum of that well. The same is 
true for p waves. However when s- waves are incident on a barrier this occurs for critical potentials which correspond to 
the supercritical values for p- waves in the spherical well (i.e. when the p-wave is at E = —m). Similarly scattering p- 
waves are in resonance when the s-wave is supercritical. This phenomenon is known as crossing and was first discussed 
by two of us [4] in a discussion of positron scattering from heavy nuclei. Under the transformation E — > — E, V — > — V, 
X — ► — X; / — * <7> the coupled equations (13) are seen to be invariant. Consequently a supercritical 2pi state for a 
well corresponds to a resonant lsi state at E = m on a barrier. 

In order to calculate the lsi phase shift for scattering from a barrier wc must use this result to calculate the 
supercritical condition for 2pi particles. By strengthening the supercritical well slightly by 1%, which, through 
crossing, increases the height of the barrier, the Is i phase shift and partial cross-section are seen to be: 



So 




k[m] 




800 850 900 950 1000 



E[keV] 



FIG. 8. The Is i phase shift and energy dependence of sin 2 So for the spherical barrier (a = 8 fm and V = 79.802 MeV) . 

We see a very clear s-wave resonance with width Tq = 43.1 keV. Similarly if the supercritical well for the lsi state 
is strengthened by 1% the phase shift for scattering from the barrier V — 78.777 MeV is: 



0.5 1 l/S 2 2T5 



k[m] 



sin 2 81 



0.5 



850 900 950 1000 



E[keV] 



FIG. 9. The 2pi phase shift and energy dependence of sin 2 5\ for the spherical barrier (a — 8 fm and V = 78.777 MeV). 

The width of the lsi resonance is seen to be Ti = 11.8 keV and the ratio of widths of the Is 1 -wave to the 2pi-wave 

2 2 2 

is found to be 3.7. 



12 



GENERAL POTENTIALS 



The reasoning we have used to understand the analytical results obtained for the spherical well applies equally to 
monotone short-range potentials with more general shapes that vanish at infinity. We shall suppose that the potential 
V(r) is spherically symmetric and is given by 

V(r) = ± vw(r/a) (68) 

where v > is a depth parameter, a is a range parameter, x — r/a > is a dimcnsionless distance, and w(x) is a 
dimensionless shape function. Thus w(x) = e~ x represents a Gaussian potential and w(x) = 1, x < 1; w(x) = 
0, x > 1 is the spherical well which we have already studied analytically. Clearly we cannot always expect to find exact 
analytical solutions. Numerical methods are flexible and powerful but involve arbitrary features, such as iteration 
method, step size, and boundary at 'infinity', which introduce corresponding uncertainties. Our approach is therefore 
to build a programme, a kind of 'virtual laboratory', to test it on the spherical well, and then to apply it to the 
Gaussian potential. For V(r) < 0, this programme would have to be able to find discrete eigenvalues; for potentials of 
either sign, it would need to be able to compute the wave-function amplitudes at the origin resulting from spherical 
waves of unit amplitude which arc directed towards the origin from far away. Since the scattering process is, in 
principle, reversible, we shall actually start waves with amplitude C at the origin, and then adjust C so that the 
waves at distant points have unit amplitude. In addition, we shall need to define and compute a 'numerical phase 
(shift)'. This question has been discussed in detail in reference [13]. 

We shall now explain more fully what we mean in this discussion by the terms 'amplitude' and 'numerical phase'. 
For very small values of r, the potentials we study are essentially indistinguishable from spherical wells or barriers. 
Hence the boundary conditions at the origin are the same as those for the exact spherical well solutions (Eq. 14). We 
shall find it convenient to separate the sign and the magnitude of \ by use also of the symbols k = \\\ = j + |, and 
T = x/lxl = x/k- In the small-r region we may write for x > 

/ = Cr k+1 

g = C'r k = aCr k 
rg _ 2k + 1 _C 
T ~ E + m-V ~ ~C 



= a 



y/p + g 2 = \C\r k Vr 2 + a 2 » a\C\r k , r<l (69) 

and for x < 

/ = Cr k 

g = C'r k+1 = aCr k+1 
g _V-E + m_C'_ 
rf ~ 2k + 1 ~ ~C ~° 
VP+9 2 = \C\r k ^l + r 2 a 2 w \C\r k , r < 1. (70) 

Ignoring the /g-ratio factor a, we shall say of these cases that the amplitude at the origin is C. We therefore seek the 
value of C which yields \J p + g 2 = 1 as r — > oo. We recall that the 'radial functions' f(r) = rip(r) and g(r) = r<j>{r) 
each have in their definitions a factor r; for bound states we have (omitting the angular factor Air) the corresponding 
normalisations, 

/ y/P{r)+g*(r)dr = / yj i\) 2 (r) + <p (r)r 2 dr = 1. (71) 
Jo Ja 



If we consider the case j = 1/2 =4- k = 1, then we have for x > y / ^' 2 (0) + <p 2 {Q) = c\C\, and for x < 
\/V 2 (0) + </> 2 (0) = |C|. This example explains our use of the phrase 'amplitude at the origin' to describe C. 

At large distances the asymptotic scattering solutions are essentially sinusoidal, like those of the spherical barrier 
(Eqs. 26,27). In this region, we choose a certain node v in g, thus g(r v ) — 0, and the amplitude \fp + g 2 at this 
point is given by |/(r„)|. Since for the Gaussian potential, the potential value V(r) is never zero, we must choose a 
sufficiently far node that the 'sinusoidal' asymptotic region has essentially been reached. For the spherical well and 
Gaussian potential we have found v = 20 to be satisfactory. Our task then is to fix v and then determine C so that 
| /(»"„) | = 1. Equivalently, we can set C = 1, and then use the resulting \f{r l/ )\~ 1 as the value of C sought. It may be 
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helpful at this point to see the wave functions in a particular case: they are shown in Figure 10 for the scattering of 
si -state waves of momentum p = 0.1 by a Gaussian barrier. 




(a) The large (/) and small (g) components for r < 7. (b) The 'large component' / has a maximum of 1 near the 

20th node r2o = 635.2 of the 'small component' g. 

FIG. 10. Resonant s-wave wave function components for a particle with incident momentum p = 0.1 scattering from a 
Gaussian barrier V(r) = v exp(— r 2 ), v = 6.8. 

The node r2o of g has been detected at r2o = 635.2, and the corresponding maximum of / has the value 1 with the 
choice C = 12.8. The final determination of the maximum is achieved by fitting a sine function to /(r) locally, as will 
now be explained. 

The position of the maximum allows us at the same time to determine what we shall call the 'numerical phase' 6. 
We use a modular function to compare the position of the peak with respect to a grid. The details of our function have 
some arbitrary elements which have been designed to yield a suitably scaled indicator of the phase. We start with 
an explicit model of the wave function for the large-r asymptotic region. The step size of our numerical integration 
procedure is represented by h. The i/ th node r = r„ of g indicates the approximate position of the maximum in /. We 
write 

f(r)=a = Dsm(pr + S'), = f(r + h), 6 = pr + 5' . (72) 
Provided a + 1 (3 and sin(#) ^ (special equations are required for the rare cases of equality), we find 

, / sin(p/i) \ a 

9 = t£m ( h\ Rl ' D = ~W\ ■ 73 

\cos[ph) — pi a J sm(fl) 

We now define 5 by the equations 

S 1 = e-pr + (l-r)^, 8 = 5 1 -w[6 1 /n\ + ^, (74) 

where we have used [^J to represent the integer part of x (often written as int(x) in computer languages). This 
definition of 5 yields a number in the range 5 € [—§,§]■ We have added a plot of S(v) to the top of our resonance 
graphs of C(v), presented below. Eq. (73) also indicates explicitly how C(v) is found: the value of v is chosen; the 
value of C in Eq. (70) is set to 1; the value of D is determined by Eq. (73); then we have C(v) = D^ 1 . In order to 
plot Figure 2b, the value of C in Eq. (70) is now set to C(v) just found, and this leads to D = 1, as the graph shows. 
All this computational activity is packaged in the form of a C++ class called 'dirac'. An instance 'dc' of dirac is 
created and messages are sent to dc to set the potential and the coupling v and the initial conditions etc; dc is then 
requested to compute eigenvalues and to perform scattering 'experiments', and to report the results. 
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This scattering quantity can be written in a closed form for the spherical potential. We start by looking at the 
asymptotic form of the wave function for very large values of r given by Eqs. (26 to 30), it is possible to find a value 
of r = R where the wave function has a peak (or node) in the top component and a corresponding node (or peak 
peak) in the bottom component. One such set of values is given by 



where N is an arbitrarily large integer. From this choice the / component will exhibit a node whilst the g component 
will be at a peak or trough of the sinusoidal wave. The modulus of the wave function at r = R from Eqs. (26 and 
27) is seen to be 



1 y u E+m E+m K ' 

p is a measure of the amplitude of the wave function at the origin when the modulus of the incident wave at great 
distances is unity. It is therefore important to examine the form of the wave function very close to the origin. The 
components of the wave function inside the potential have the form given by Eqs. (14). Using the following identity 
for the limiting value of the spherical Bessel function for small arguments 

lim [x- n j n (x)] = J^~frj^. (7?) 



*(r) 



f(r) 



ai 



r(pr)'x 

(2i x + l)!! 

(pr)'-v+i 



(78) 



J^E-V+m (2Z_ X + 1)<! 

The modulus of the wave function at the origin is therefore 



«i 



P 



x 



l*(0)| = VFW+gW) 



E- V + m (2x- 1)!! 

P-*- 1 ..-V 



ai 



(-2 X -1)!! 



X>0 
X<0 



(79) 



In particular for s-waves 



whilst for p-waves 



|*(0)| ~ air 



|*(0)| 



ai 



pr 



E-V 



(80) 



(81) 



in agreement with Eqs. (69 and 70). Following the argument given above we can set a\ 
where the sign indicates a barrier (+) or well (-) - is found by 



C(±) = 



1 and the quantity C(±) 



(82) 



This can be established upon substitution of Eqs. (76, 23 and 24) into the above. C(— ) can be established from the 
above by the transformation V — > —V. 

For vn — 1 we decided first to study three resonance curves, that is to say, graphs of C(v) for a range of v > 
large enough to reach just beyond the third resonance (this is illustrated for the Gaussian potential). We kept j = \ 
constant and studied the four cases: r = — 1 (si), r = 1 (pi), repulsive barrier (+), and attractive well (— ). For both 
of the potentials we found that the resonances associated with the attractive wells corresponded exactly to the critical 
eigenvalues at E — m in the same sector, s or p; whereas, for scattering from barriers, the resonances corresponded, 
in agreement with our crossing theorem, to the eigenvalues E = — m in the 'other' sector, s to p, and p to s. For the 
spherical potential, the exact critical couplings are given by the following formulas in the p-sector: 
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pi(-) => v = (1 + (rnr) 2 ) 2 - 1, n= 1,2,3,.. 
pi(+) w=(l + (n7r) 2 )'+l, n= 1,2,3,.. 



(83) 
(84) 



which can be established from Eqs. (47 and 51) with a = 1. The corresponding values in the s-sector are found as 
solutions (Eqs. 46 and 53) to a quasi-bound state problem in which the spinor is not normalised since <j>(r) approaches 
a constant non-zero value as r — > oo. The numerical values relevant to the resonance curves are shown in Table 1. 



9. V ' 


si(-) 




fiH 


5.27 
8.40 
11.54 


1.11 
4.20 
7.33 


4.30 
7.36 
10.48 


2.30 
5.36 
8.48 



Table 2: Critical couplings for scattering from spherical barriers (+) and wells (— ). 

As we have explained, these critical couplings correspond to (V < 0) eigenvalues E = ±m by the rules: 

~(E= -m)pi, 
pi(+) ^(E= -m)si, 
si(-) ^(E = m)si, 
pi(-) = m)pi. 

The corresponding C(p) resonance curves are shown in Figure 11. 



(85) 
(86) 
(87) 
(88) 




(a) C(p) for s-wave scattering from a spher- 
ical barrier of height U = 5.29. 



(b) C'(p) for p-wave scattering from a spher- 
ical barrier of height U = 4.4. 
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(c) C(p) for s-wave scattering from a spher- (d) C'(p) for p-wave scattering from a spher- 

ical well of depth U = -4.195. ical well of depth U = -2.2. 

FIG. 11. These plots show the scattering measurement C(p) plotted against the incident momentum p for s and p waves 
incident on supercritical spherical barriers and subcritical wells. 

By repeating the analysis for the Gaussian potential, we obtain the resonance peaks and corresponding eigenvalues 
E = ± m at the critical couplings shown in Table 2. 



9. V ' 


si(-) 




PiH 


6.75 
10.42 
14.04 


1.26 
4.37 
7.70 


5.62 
9.23 
12.83 


2.96 
6.11 
9.44 



Table 3: Critical couplings for Gaussian barriers (+) and wells (— ). 

The resonance curves C(v) for an incident momentum of p = 0.1 for the Gaussian potential are shown in Figure 12. 
We have also computed the critical couplings corresponding to s and p states with E = —m = — 1 and E = in = 1. 
Allowing for s <-> p crossing in the E = —m cases, these values agree with the scattering results reported in Table 1. 
The heights of the peaks in these resonance curves are somewhat arbitrary: they depend on the choice of p, and, more 
randomly, on the number of v values used in the plot. Ideally we should have p — > : the barrier peaks are sharper for 
smaller p, and their heights become unbounded; these effects can be 'seen' numerically by fine explorations of C(v) for 
v near critical values. However, as p is reduced, the scale of the relevant part of the wave function expands, and the 
numerical integration must be carried out to greater and greater distances. The choice p — 0.1 and its implications 
are therefore arbitrary. Moreover, the plots we exhibit were made with n v = 400 points in the abscissae: the closeness 
to which the resulting values of v lie to critical values is therefore a matter of chance. The resulting numerical maxima 
have these arbitrary aspects, which, of course, emphasise, by comparison, the great importance of the study of exact 
analytical solutions when such are available. 
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(a) s-wave resonances and phase shift for the 
Gaussian barrier. 



(b) p-wave resonances and phase shift for 
the Gaussian barrier. 




J 


I 

.1. J 





(c) s-wave resonances and phase shift for the 
Gaussian well. 



(d) p-wave resonances and phase shift for 
the Gaussian well. 



FIG. 12. These plots show the scattering measurement C(p) plotted against the incident momentum p for s and p waves 
incident on supercritical Gaussian barriers and subcritical wells. 

We have found similar results for other short-range potential shapes, such as the exponential potential w(x) — e~ x 
and the Wood-Saxon potential w(x) = 1/(1 + e x_1 ). 



CONCLUSIONS 



We have reviewed the concept of a resonance in potential scattering in quantum mechanics in terms of the Wigner 
criterion that not only must the phase shift pass through n/2 but there must also be a positive time delay in the 
scattering process. In non-relativistic scattering from a monotone spherically-symmetric potential well that implies 
that s-wave resonances do not exist, but that p-wave and higher angular momentum resonances do exist. The zero 
energy solutions of the Schrodinger equation provide a starting point for the study of resonances for scattering off a 
well: resonances are obtained by slightly weakening the potential needed to obtain a zero energy solution. 

In the case of relativistic potential scattering using the Dirac equation, the zero energy solutions are replaced by 
zero momentum or critical solutions at E — m and E = —m. Now resonances are obtained by slightly weakening 
the critical potential well giving E = m or slightly strengthening the supercritical well giving E = —to. In the 
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case where a critical solution at E = m exists for a potential well, slightly weakening the potential will give rise to 
resonances in s-waves as well as in higher waves. Slightly increasing the strength of the the supercritical well will 
give rise to resonances in scattering from a potential barrier using the crossing theorem, which relates positive energy 
solutions of the Dirac equation for a potential barrier to negative energy solutions for a potential well. We thus find 
that resonances exist for scattering from barriers in s-waves and higher waves. In particular crossing relates s-wave 
solutions for a well to p-wave solutions for a barrier and vice versa. 
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